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ABSTRACT 

We present high spatial resolution (< O'.'S ~ 40 AU) Submillimeter Array observations of the 865 fim 
continuum emission from the circumstellar disk around the young star DoAr 25. Despite its bright 
millimeter emission, this source exhibits only a comparatively small infrared excess and low accretion 
rate, suggesting that the material and structural properties of the inner disk may be in an advanced 
state of evolution. A simple model of the physical conditions in the disk is derived from the sub- 
millimeter visibilities and the complete spectral energy distribution using a Monte Carlo radiative 
transfer code. For the standard assumption of a homogeneous grain size distribution at all disk radii, 
the results indicate a shallow surface density profile, S c>c r~P with p « 0.34, significantly less steep 
than a steady-state accretion disk {p = 1) or the often adopted minimum mass solar nebula {p — 1.5). 
Even though the total mass of material is large (M^ « 0.10 M0), the densities inferred in the inner 
disk for such a model may be too low to facilitate any mode of planet formation. However, alternative 
models with steeper density gradients (p « 1) can explain the observations equally well if substantial 
grain growth in the planet formation region (r < 40 AU) has occurred. We discuss these data in the 
context of such models with dust properties that vary with radius and highlight their implications for 
understanding disk evolution and the early stages of planet formation. 
Subject headings: 



1. INTRODUCTION 

The details of the planet formation process in its early 
stages are largely determined by the physical conditions 
in the disk around the host star. The densities, temper- 
atures, material content, and structure of this gas and 
dust reservoir critically influence the mechanisms and ef- 
ficiencies of the assembly of a planetesimal population, 
as well as its subsequent growth and dynamical evolu- 
tion. Spatially resolved observations of disks at millime- 
ter wavelengths provide unique access to these physical 
conditions, inclu ding in particular th e spatial distribu- 
tion of mass fe.g-. IWilner fc Lavll2000 f). Such constraints 
on the physical structure of disks around young (^1 Myr) 
T Tauri stars are expected to offer a glimpse at the rela- 
tively p ristine initial conditions av ailable for planet for- 
mation (jBeckwith fc Sargentlll996l ). 

However, recent observations of some disks raise ques- 
tions about how pristine or initial these conditions really 
are, even at such young ages. A small collection of "tran- 
sition" disks show pronounced dips in their infrared spec- 
tral energy distributions (SED) which sugg est that large 
inner hole s or gaps have b een cleared (e.g.. ICalvet et al.l 
[2005; Esp aillat et al.ll2007|)i a scenario confirmed in some 
cases by direct imaging (jPietu et al.ll2006t iHughes et al.l 
l2007HBrown et al.ll2008l ). Although less dramatic, a sep- 
arate population of disks with relatively faint infrared 
SED excesses (although without distinct SED gaps) but 
bright millimeter emission can be interpreted as evidence 
for substantial particle growth — an d thus diminished 
grain emissivity — in the inner disk (|Lada et all l2006l : 
iNaiita et al.ll2007b . These indicators suggest that some 
fraction of ~1 Myr-old disks have already made signifi- 
cant progress along the path toward making planets. 

In this Letter, we use high spatial resolution submil- 
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limeter observations to help interpret the physical condi- 
tions in the disk around DoAr 25. Loca ted in the L168 8 
dark cloud in Ophiuchus {d sa 145 pc; lMakarovll2007f) . 
the DoAr 25 disk has a comparatively small infrared 
excess and accretion rate (M ~10"^° -10^^ Mf?) yr~^; 
iMuzerolle e t al.l ll998HNatta et al.ll2006l) but b right mil- 
limeter emission (|Andrews fc Wilhamsl r2007b[ ) . Follow- 
ing a description of the observations and their calibration 
in §2, we utilize a standard disk model and radiative 
transfer calculations to estimate the disk structure and 
mass distribution in §3. We discuss the results in the 
contexts of particle growth and planet formation in §4. 

2. OBSERVATIONS AND DATA REDUCTION 

DoAr 25 {a = 16'^26'"23!68, 6 = -24°43'14'.'l [J2000]) 
was observed in the "very e xtended" config uration of the 
Submillimeter Array fSMA jHo et al.ll2004 fl on 2007 May 
26 and June 17 in exceptional conditions, with zenith op- 
tical depths of <0.06 at 225 GHz. In this arrangement, 
the 8 SMA antennas (6 m diameter each) span physical 
baselines up to 508 m. The observations cycled between 
DoAr 25 and two quasars (J1517-243 and J1626-298) 
on 10 minute intervals. The raw visibilities were cali- 
brated with the MIR package. Passband calibration was 
conducted with bright quasars (3C273, 3C279), and com- 
plex gain calibration was performed with J1626— 298. 
MWC 349 and Callisto were used to set the absolute 
flux scale, which is accurate to ~10%. Data from both 
nights and each of the two 2 GHz sidebands were com- 
bined. Using J1517— 243 as a consistency check on the 
phase transfer, we estimate that phase noise and baseline 
errors generate an effective seeing of <0'.'l. 

Additional observations in the compact and extended 
SMA configurations (baselines of ~15-200m) w ere ob- 
tained in 2006 (see [Andrews fc Winiami l2007a[ ). The 
data show excellent consistency between all configura- 
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Fig. 1. — (left) High resolution SMA image of the 865 fim continuum emission from the DoAr 25 disk. Contours start at 7.5 mjy beam ^ 
(3cr) and increase by factors of \/2. The synthesized beam size is shown in the lower left corner, (middle) The elliptically-averaged 865 /^m 
visibilities (the real part) as a function of the deprojected baseline length, (right) The DoAr 25 SED, where the ordinate corresponds to 
4TTd?uF^. Optical photometry is from lVrba eTall {TgH, UBV) and lWilking et al.l ||2005|. RI). Inf rared data are from 2M ASS ( Cut ri et all 
120031) . the Spit zer Space Telescop e (L. E. Allen 200 7, private com munication) , 'BaT:son Y et al.l ||2005I1 . and lBontemps et al.l ll200lV Millimeter 
fluxes are from lDent etaLl 1 11993 ). this paper, and [Andre &: Montmcrlc (1994). Error bars include uncertainties in the absolute flux scales. 
The grey curve shows the Kurucz model adopted as the input stellar radiation source for the radiative transfer modeling. Overlaid on the 
visibility and SED panels are the model disk behaviors for a standard grain size distribution (red: §3) and an alternative distribution (blue) 
with grain growth (and reduced emissivity) in the inner disk (r < 40 AU; §4). 



tions on overlapping spatial scales. All continuum data 
were combined to provide visibilities at an effective fre- 
quency of 346.9 GHz (865 /xm). For the 2006 data, we 
simultaneously observed the CO J=3— 2 transition at 
345.796GHz with a spectral resolution of 0.70km s~^. 
The MIRIAD package was utilized for the standard 
tasks of Fourier inversion, deconvolution, and imaging 
of the calibrated visibilities. A submillimeter continuum 
map was generated with a Briggs robust=0.25 weight- 
ing scheme, providing a synthesized beam FWHM of 
0'.'43 X 0'.'32 at a position angle (PA) of 15°. 

3. RESULTS 

Figure [1^ shows a high spatial resolution 865 /im con- 
tinuum image of the DoAr 25 disk. An elliptical Gaussian 
fit to the visibilities indicates a well-resolved source with 
an integrated flux density of 563 ± 3 mJy (not including 
the calibration uncertainty) and an inclination of 62 ±3° 
at a PA of 111 ± 3°. The SMA visibilities are presented 
in Figure [TJd, cUiptically averaged afte r their deprojec - 
tion based on the above geometry (see iLav et al.lll997D . 
The visibilities drop off steeply, suggesting a shallow gra- 
dient in the radial surface brightness profile of the disk. 
The SED for DoAr 25 is shown in Figure [T];, using data 
compiled from the literature. A relatively small thermal 
excess is noted above the stellar photosphere out to mid- 
infrared wavelengths, before the SED picks up to show 
some of the brightest mi Uimeter emission of all the T 
Tauri disks in Ophiuchus () Andrews fc Wilham"3l2007b f). 

For a quantitative view of the physical conditions in the 
DoAr 25 disk, we attempted to reproduce the full SED 
and SMA visibilities with the 2-D axisymmetric Monte 
Carlo radiative transfer code RADMC (v3.1, C. P. DuUe- 
mond 2007, private comm unication). This code utili zes 
the algorithm developed bv lBiorkman k WoodI ((2OOTI to 
compute temperatures in a circumstellar medium heated 
solely by stellar irradiation and with a density structure 
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where the radial coordinate (r) runs from the dust sub- 
limation radius to an outer boundary (Rd), the angular 
coordinate (9) runs from the rotation axis {9 — 0) to the 



midplane {9 — 7r/2), and the surface density and scale 
height vary with radius as power laws (S — Y,o{r / Rd)~^ 
and H = Ho{r / Rd)'^)- Aside from some minor coor- 
dinate conventions, this disk structure model is iden- 
tical to others commonly used in the literature (e.g., 
[Chiang fc Goldreichlll997|). The code itself was described 
by DuUcmond & DominiH (|2004f) . but now includes a dif- 
fusion algorithm that is employed for grid locations with 
poor photon statistics (i.e., near the disk midplane). 

The central illumination source was a Kurucz model 
(gray curve in Fig. Ic) fixed to match the proper- 
ties of the star and optical SED: a K5 spectral type 
(Teff= 4250K), L, = 1.3Lo, and Av^2.9 (IWilking et all 
|2005[ ) . We initially adopted the emissivity properties 
provided bv lD'Alessio et all (|2001h . for a dust grain size 
(a) distribution that varies as n(a) oc a~^'^ between 
amin=0.005/.tm and amax=lmm. These emissivities are 
essentially identical to the standard disk values used at 
millimeter wavelengths (J3cckwitii_ ct ah 1 990t ) . With this 
information, we ran the radiative transfer code over a 
grid of 5 free parameters that characterize the disk struc- 
ture: Rd, p, h, and the normalizations of S (i.e., the disk 
mass, Md) and H. For each point in the grid a model 
SED and 865 /im visibilities were generated with a ray- 
tracing program, assuming the above orientation. 

The best-fit model parameter values estimated from 
minimizing are Rd = 310 AU, p — 0.34, h — 1.11, 
Md = O.IOM0, and H = 28AU (at Rd). The synthetic 
SED and visibilities for a model with these parameters 
are shown in red in Figure 1. The reduced value com- 
paring the best-fit model and the data shown in Fig. lb 
and Ic was ~2, primarily due to difficulty reproducing 
the mid-infrared spectrum in detail. The formal 1 a pa- 
rameter uncertainties are '--^15AU for Rd, '--^O.OIMq for 
Md, and ~5AU for H (at Rd), while p can range from 
0.2—0.4 and the flared shape of the disk can run as low as 
/i « 1.07 (but not higher than the best-fit value without 
over-producing infrared emission) . 

As an a posteriori test of the derived disk properties, 
Figure[2]compares an optical (0.6 /im) Hubble Space Tele- 
scope {HST; WFPC2/F606W) image of starlight scat- 
tered off the DoAr 25 disk surface (jStapelfeldt et al.l 
l2008f ) with a RADMC model prediction corresponding to 
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Fig. 2.— (ie/t) An optical HST image (F606W, ~0.6^tm) of 
starlight scattered ofT the surface of the DoAr 25 disk, (right) The 
scattered light image predicted from the best-fit model parameters 
described in §3, based solely on the SMA visibilities and the SED. 
The stellar position is labeled with a star symbol, a gray line marks 
the orientation of the disk major axis, and a logarithmic bright- 
ness scale is shown on the right. While the model does not include 
direction-dependent scattering (and is therefore not used in the 
structure determination), the comparison with the HST data con- 
firms the size scales, orientation, and general disk geometry that 
were derived from the SMA data alone. 

the above parameters. Because this version of RADMC 
only treats isotropic scattering, a quantitative compari- 
son with the data is unjustified. However, it is clear that 
the best-fit model is in good qualitative agreement with 
the disk orientation and geometry as traced by both the 
scattered flux scale and the location and intensity of the 
dust lane marking the midplane. 

While not shown here, channel maps of the CO J— 3— 2 
transition from the SMA data (with a I'.'S x f'.'6 beam 
at PA = 4°) exhibit two distinct patches of weak line 
emission separated by ~4km s~^ toward DoAr 25, cen- 
tered on the local cloud velocity {Vlsr = 3-3 km s~^). 
Single-dish mapping of the Ophiuchus clouds in the CO 
J=l— transition shows a bright, optically thick line 
with a FWH M of ^3 km s~^ coincident with the DoAr 
25 position ( Ridge et al.ll200^ . Unfortunately, this sug- 
gests that the line emission from the disk remains sub- 
stantially contaminated by confusion with the more ex- 
tended molecular cloud. 

4. DISCUSSION 

The spatial distribution of mass, characterized by the 
surface density index p, is a critical diagnostic for un- 
derstanding the material evolution and planet formation 
potential of a disk. It reflects the disk viscosity and there- 
fore provides an observable indication of the evolution 
timescale dictated by the accretion process and angular 
momentum conservation. Moreover, the mass distribu- 
tion will influence the composition and architecture of a 
developing planetary system. By comparing the decay of 
accretion rates over time with models t hat use a simple 
prescription for the turbulent viscosity, iHartmann et al.l 
argue for a typical value p=l. A steeper gradient 
(p=1.5) is invoked to account for the mass distribution in 
the solar system, once augmented to cosmic abundances 
and smeared into^ annuli (t he minimu m mass solar neb- 
ula, or MMSN: lWeidenschillingjll977[ ). 

The low value of p derived for the DoAr 25 disk would 
have impor tant implications for un derstanding its evolu- 
tion. In the lHartmann et al.l (|1998f ) prescription for a vis- 
cous accretion disk, such a low value of p would produce a 
significantly diminished mass accretion rate onto the star 
over most of the disk lifetime, perhaps lower by a factor 
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Fig. 3. — (top) The radial surface density profile (S; gas+dust) 
for the best- fit model described in §3 (red) and the toy model with 
large grains in the inner disk described in §4 (blue). The latter is 
decomposed into two grain population s, with am ax = 1 mm outside 
40 AU and 1 m (both from D'Alessio e t al.l2001l ) inside that radius. 
The minimum mass solar nebula (MMSN) density profile is shown 
in black for comparison, (bottom) The radial profile of the 865 /im 
continuum optical depth (t = ftS; where k represents the material 
emissivity) for the same models. The open diamonds track the 
optical depth that would be expected if the standard (amax = 
1 mm) grain size distribution emissivity were adopted everywhere 
in the disk for the steep density gradient (p = 1) model. The SMA 
resolution limit for these data is marked with a vertical dotted line. 

of ^2 compared to the canonical disk with p — 1. This 
is in reasonable agreement with c onstraints on M, which 
vary from low (5 x 10"^ Mq vr- MMuzerolle et al.lll998D 
to ne gligible values (< 5 x 10~^° M^ vr~^: iNatta et all 
[20061) . In essence, the accretion process is less efficient 
in a disk with a shallow density gradient. The same 
could be said about the formation of planets in such a 
disk. In spite of the large M^, the inner disk densities 
implied by this low-p model are more than an order of 
magnitude smaller than those required b y planet forma- 
tion models (e.g.. iHubickvi et al.ll2005l also note that 
even though M^/M* « 0.18, the Toomre stability cri- 
terion holds at all radii in this model). Previous mil- 
limeter observations of disks have general ly found large 
p (> 7; e.g., ;W ilncr ct al. 19 91 [Andrew s & Wi lliamsl 
I2007at iPietu et al.i i200& iHamidouche et al.i i2006l ). al- 
though the derived value for DoAr 25 lies in the low end 
of a wide distribution; the variety of modeling techniques 
and spatial resolutions in those studies make direct com- 
parisons difficult. 

Of course, there are alternative interpretations of the 
observed DoAr 25 disk properties that do not call for such 
a small density gradient. The absence of a pronounced 
dip feature in the infrared SED or an apparent "hole" in 
the submillimeter surface brightness distribution can ef- 
fectively rule out a substantially cleared inner region like 
those those in transition disks. However, because the 
millimeter continuum emission traces the optical depth 
in the disk — the product of the density structure and the 
dust emissivity — the observed shallow brightness profile 
could be explained equally well with steep density gradi- 
ents {p > 1) if the emissivity compensates by decreasing 
toward small radii. For the models described in §3 (and 
generally in the literature), the grain size distribution is 
assumed to be spatially uniform throughout the disk. In 
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reality, particle sizes in the inner disk might be larger 
than in the outer disk because the coUisional growth 
timescale is shorter for high er densities and orbital ve- 
locities (^e.g.. iDullemond fc D ominik 2005; Garaud 2007, 
and references therein). ^'Alcssio et al.l (|2001l) show 
that the millimeter emissivity quickly decreases as amax 
grows beyond '^l mm. Therefore, a shorter grain growth 
timescale in the inner disk would naturally produce an 
emissivity that decreases at smaller radii. 

To illustrate the effects this would have observation- 
ally, we constructed a crude toy model that has the stan- 
dard steady-state accretion disk density profile {p = 1), 
but allows particles with a factor of ~5 smaller emis- 
sivity (at 865 /xm; corresponding to Cmax ~ 1 m) to ac- 
count for essentially all of the mass on solar system scales 
(r < 40 AU). As shown in blue in Figure [U models that 
mimic grain growth can reproduce both the DoAr 25 
SED behavior and the shallow millimeter brightness pro- 
file fairly well. However, given the additional freedom 
that comes with more model parameters, the data do not 
uniquely distinguish values of inputs like Omax- Differ- 
ent values can be accomodated by adjusting the relative 
abundances and locations of the grain populations. Fig- 
ure [3] shows the radial distributions of S and the 865 fiui 
optical depth for both this toy model and the more stan- 
dard prescription adopted in §3. While the details in 
such an ad hoc demonstration model are unimportant, 
the lower panel of Figure [3] highlights the salient point: a 
decrease in the inner disk emissivity can effectively mask 
a steep density gradient. 

There is considerable evidence from the shape of mil- 
limeter continuum spectra that disk emissivities have 
been modified by substantial particle growth, albeit 
in an average sense due to limited spatial resolution 
(e.g.,LB eckwith fc Sargent"1991'; 'Testi et al."2001|, [20031 : 
IWilner et al. 2005; Rodmann ct al. 2006). DoAr 25 has a 
fairly typical long- wavelength SED shape (n — 2.7 ±0.2, 
where F^, cx i^") for a T Tauri star. High spatial res- 



olution observations at multiple frequencies (from ~40- 
400 GHz) can potentially disentangle the millimeter SED 
shape (n; i.e., colors) in the evolved inner disk from the 
more primordial material at larger radii. Such data could 
be used to empirically reconstruct the shape of the emis- 
sivity spectrum as a function of location in the disk, 
which could subsequently be modeled to map out the 
constituent grain populations. 

If the properties of disks like DoAr 25 are caused by 
particle growth in the inner disk, they can be considered 
less-evolved versions of "transition" disks, with larger 
remnant optical depths in their inner holes (effectively 
muting out the signature dip in the SED noted for tran- 
sition disks). It would be interesting to model high- 
resolution millimeter observations of similar disks and 
see if shallow density gradients (low p) are commonly in- 
ferred when assuming a standard grain size distribution. 
If that is the case, such results could be a signature that 
a few percent of ~lMyr-old disks are already well on 
their way to making planets. 
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